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Abstract 

We discuss exclusive electromagnetic production of two neutral p mesons in coherent photon- 
photon processes in ultrarelativistic heavy-ion collisions. The cross section is calculated in the 
equivalent photon approximation (EPA) . Both uncertainties of the flux factors and photon-photon 
cross sections are discussed in details. We show that inclusion of precise charge densities in nuclei 
is essential for realistic evaluations of the nuclear photon-photon cross sections. We find that 
the cross section, especially with realistic flux factors, is sensitive to low energy in the subsystem 
77 — ► p°p°. The experimental data for the 77 — ► p°p° cross section extracted from e + e~ collisions 
are parametrized and used to estimate the nucleus-nucleus cross section. In addition, we include 
vector-dominance-model(VDM)-Regge contribution which becomes important at large photon- 
photon energy. Large nuclear cross sections are obtained. We discuss a possibility of focusing 
on the large-energy component. We find that both p° mesons are produced predominantly at 
midrapidities and could be measured by the STAR collaboration at RHIC. 
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I. INTRODUCTION 



Exclusive production of elementary particles (lepton pairs, Higgs, etc.) or mesons (vector 
mesons, pair of pseudoscalar mesons, etc.) in ultrarelativistic heavy ion collisions is an 
interesting and quickly growing field 0, 0, 1 of theoretical investigation. On experimental 
side the situation is slightly different. So far only single-p° exclusive cross section AA — > 
AAp° was measured [If. Here the dominant mechanism is the photoproduction described 
by the photon-pomeron fusion. 




FIG. 1: The reaction discussed in this paper. 

In the present paper we consider exclusive production of p°p° pairs in 77 collisions. So far 
only integrated cross section for A A — > AAp°p° was roughly estimated in the literature 
using a pQCD-inspired modelling of the Pomeron exchange between mesons for AA — > 
AAJ/tyJ/ty and A A — > AAp°J/ty and assuming Regge factorization. Uncertainties of 
almost three-orders of magnitude were found. As we discuss here, in the AA — > AAp°p° 
reaction one is sensitive rather to the close-to-threshold region for the 77 — > p°p° subprocess 
where the specific pQCD approaches of [HJ are not applicable at all. Furthermore the Regge- 
factorization for this hard and soft processes is likely broken. The p°p° pair will eventually 
decay into the final state of four charged pions. If the p°p° invariant mass is very large, the 
two 7r + 7r~-pairs will be separated by a large rapidity gap. It has been recently suggested 
that in this kinematics, charge asymmetries can be used to discover /study the QCD-Odderon 
exchange 

At the RHIC energies discussed here, we will be however dominated by 77 center-of-mass 
energies in the region of a few GeV. Here, experimental data for 77 — > p°p° are known 
and were measured by several groups at e + e~ colliders 0. Experimentally one observes 
a huge enhancement close to the threshold. The origin of this enhancement was never 
understood. Some speculations are discussed in Ref.[8(. The huge cross section at W ~ 
1.5-2 GeV was interpreted as a tensor resonance decaying into p°p° channel j^j. In Ref. 13 



the Regge factorization was tested. In the following we shall see the consequences of the 
anomalous behaviour seen in the e + e~ data for the two-photon production of the p°p° pair 
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in the nucleus-nucleus collisions. 



II. CALCULATION 



A. Equivalent photon approximation for nuclei 



Let us consider the process AA — > AAp°p° depicted in FigfTJ The cross section takes the 
familiar form of a convolution of equivalent photon fluxes and 77-cross sections: 



da{AA -> p°p°AA; s AA ) 

— = dn ry {x 1 ,X2, 0) (7(77 -> p p ; x^saa) + ... (2.1) 

Here we omitted helicity-dependent pieces, see below for more details. The effective photon 
flux is expressed through the electric field strengths of the ions and reads (see also [l2| ) : 



dn^{x x ,x 2 ,b) = ( d 2 b l d 2 b 2 SlM5( 2 \b-b 1 + b 2 )^^\\E(x 1 M)\ 2 \E(x 2 M)\\ 

J X\ X 2 7T 

(2.2) 

where the electromagnetic field strengths are given in terms of the charge form factor of the 
nucleus F ern ((f) by: 



E(x, b) = Z^4ira em J 



d 2 q 



-ibq] 



(2tt) 2 q 2 + x 2 M\ 



expi 

dq 2 {bq)Mbq) 
4tt t ? 2 + x 2 M| X "" 



[q 2 + x 2 M\) 



F em (q 2 + x 2 M 2 A ) . 



(2.3) 



In practice F em ((f) is obtained as the Fourier transform of the nuclear charge densities 
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extracted from electron-nucleus scattering data, for a useful compilation, see [13|, 
putting F* m = 1 one obtains the well-known result for a point-like charge (see e.g. the 
textbook [151]): 



E pt (x, b) = -i Z ^ a *™ * {xM A b)K x {xM A b) . (2.4) 

Notice that our formula for the effective Weizsacker -Williams flux of photons includes 
the impact-parameter dependent absorption factor S 2 bs (b). It represents the probability 
that no inelastic interaction between the nuclei occurs. At RHIC energies it could be for 
example estimated from Glauber theory using the Czyz-Maximon [l6| approximation: 



S 2 abs (b) = exp ( - off N J d 2 sT A (b - s)T A (s)) , (2.5) 

where T A (b) is the optical thickness of the nucleus. We checked, that for our purposes it 
can be well approximated by the result for a black disc: 

S 2 abs (b) = 9(b - 2R A ) . (2.6) 



3 



The presence of the absorption factor also induces an azimuthal correlation between 
impact parameters 61 and b 2 . Now, soft photons are polarized linearly in the transverse 
plane, and therefore a subtle helicity dependence arises. Strictly speaking, one should write 

da(AA - fPfPAA^ = dn \\ {xuX ^ b) & (77 _ p o p0) + d . {xu X2j b) a±{ll _ p o p0) ? 



d 2 b 



(2.7) 



where 

dn» (x 1; x 2 ,6) = / d 2 b 1 d 2 b 2 S 2 abs (b)5^ 2 \b-b 1 + b 2 )^^^-\E(x 1 ,b 1 )-E(x 2 ,b 2 )\ 2 

J X\ X 2 7T Z 

f/ ' , ' J ' / -'- ^(.r . 6 j ) < E(.r 2 .fr,)H 2 

(2.8) 



dnUx h x 2 , b) = [ d 2 hd 2 b 2 SlMd^ib -bx + b 2 )^^\\[E{ Xl , b x ) x E(x 2 , b 2 ) 

" J X-t X? 7T Z 



are the fluxes for photons with parallel (||), and perpendicular (_L) linear polarizations; a" 
and a x are the corresponding 77 cross sections. The small effects associated with this 
helicity dependence have been discussed in ll|| and will be neglected here. 



A further comment on the practical uses of the equivalent photon approximation is in 
order. It is often pointed out that the Fermi-Weizsacker-Williams idea of the equivalent 
photon flux represents indeed the archetypical parton-model concept. In this spirit, often 
flux factors of equivalent photons are calculated as for point-like particles with charge Ze, 
and the total cross section is then evaluated using a simple parton-model type formula: 



a (AA - A(p°p°)A] = [ d^du ^^ a (77 - P°p°) ■ (2.9) 

' J UJ\ UJ 2 ' ' 



When calculating the two-dimensional integral it must be checked if W 2 = Auj\uj 2 > 4m 2 . 
Here 



n[uj} 



J d 2 bN(u,b) = ~ J d 2 b\E(x,b)\ 2 ; x = (2.10) 



is the flux of photons in an ultrarelativistic charge of energy jMa- For RHIC, at the cms- 
energy a/J/vtv = 200 GeV we have 7 ~ 100. 

Notice that such a factorized form is not borne out by the general formula (12.21) which ac- 
counts for the strong absorptive corrections. In particular, one cannot meaningfully improve 
the ^-integrated Eq.f l2.9p to effectively account for strong absorption, in particular this is 
not achieved by simply restricting the impact parameter integrals f)2.10p of the individual 
photon fluxes. 

When written in terms of photon energies Ui in the nucleus-nucleus center-of-mass, for- 
mula (12.11) takes the well-known form [ll| 



a (AA - A(p°p )A; s AA ) = f d 2 b x d 2 b 2 Q (\b x - b 2 \ - 2R A ) ^^Nfa, b 1 )N(u 2 , b 2 ) 



dio\ doj 2 
U)\ io 2 

x a (77 p°p°; 4uo x u 2 ) . (2.11) 
We also use the 77 cms-energy W 77 , and the rapidity- type variable Y defined through 
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(P°,P°) 




FIG. 2: The quantities used in the impact parameter calculation. 



77 e Y 



w. 
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77 -Y 



W, 



(2.12) 
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The cr (77 — > p°p°) cross section can be calculated from models and/or taken from exper- 
imental data, if existing. Below we discuss separately the low-W 77 and high-W^ case. 



B. Low-energy 77 



p°p° 



cross section 



The cross section for this process was measured up to W 77 = 4GeV [T7J. At low energy 
one observes a huge increase of the cross section. Several possible scenarios were discussed 
in this context (see e.g. Ref. (§]). For one real and one virtual photon induced processes the 
enhancement was interpreted as due to an isotensor meson decay [? ] . We leave the difficult 
problem of the microscopic origin of the close-to-threshold bump for future studies and take 
here a pragmatic attitude of using directly experimental data. 

In FigJ3] we have collected the world data (see 0, 13] and references therein). A huge 
rise of the cross section can be seen close to the threshold. Here we shall use rather directly 
experimental data in order to evaluate the cross section in nucleus-nucleus collisions. In 
Figj3] we show our fit to the world data. One can observe a small inconsistency of the data 
measured by different groups. This means that also our parametrization has about 20 % 
accuracy. 



C. High-energy 77 — ► p p cross section 

The cross section above W 77 = 4GeV was never measured in the past. It is well known 
that the cross section for 77 — > hadrons can be well described in the VDM-Regge type 
model. Here we present a similar approach but for the specific, well defined, final state 
channel p°p°. 

In the VDM-Regge approach the amplitude for the 77 — > p°p° can be written as 

A4 77 ^ p o p o(s, t; qi, q 2 ) = C 7 ^ p oC 7 ^ p oA / lp0. p o.^ p o p o(s, t; q u q 2 ) ■ ( 2 -13) 
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FIG. 3: The elementary cross section for the 77 — * p p reaction. In the left panel we display 



the collection of the e + e~ experimental data [7], [17| and our fit. In the right panel we show our 



predictions based on the VDM-Regge model decribed in the text. For comparison we show by the 
dashed line in the left panel also the result when the form factor correcting for off-shell effect is 
ignored. 

I 2 

Above C 7 ^ p o = J P and we use f p = ySt. The later value is obtained in order to reproduce 



'y-^ff — y j p vw. uo>, Jp — 2 54 . 
the decay width of the p° meson into dileptons. 

For energies PF 77 > 2-3 GeV the amplitude for the p°*p°* — > p°p° proceses can be written 
in the Regge form : 

M p o* p o*^ p o p o(s, i) = s\ rnp(s, t) C F ( — J + rj R (s, i) C M ( — j J -F(i; qf) F(i; q\\ 

(2-14) 

Above we have introduced vertex form factors which, in general, are functions of exchanged 
pomeron/reggeon four-momentum and photon (p° meson) virtualities. We parametrize them 
in the factorized form: 

Fii-A 2 ) = ex P • exp (^^) • (2.15) 

The second term, which "corrects" for p° meson virtuality, is normalized to unity when p° 
meson is on mass shell. We expect the slope parameter of the order B ~ 4 GeV -2 and 
the parameter responsible for off-shellness of p° mesons A ~ 1 GeV. We take the powers 
of the pomeron (ap) and reggeon («jr) terms from the Donnachie-Landshoff fit to the 



1 We assume helicity conservation and neglect helicity flip in the present analysis. As far as only a good fit 
of the total cross section is concerned this is of no further relevance for our purposes. 
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total NN and ttN cross sections 18j. The parameters Cp and Cm are obtained assuming 
Regge factorization and assuming that a(p°p° — > p°p°) = cr(7r 7r — > 7r°7r°) (see e.g. 1^ , 
they are: Cjp = 8.56 mb, Cjr = 13.39 mb). While this seems justified for the pomeron 
term, it is not so obvious if it is true for the reggeon terms. Consistent with our choice 
of normalization rjp(s,t) and ?7jr(s, t) are complex functions such that: T]p(s,t = 0) « % 
and T]]fi(s,t = 0) ~ % + 1. Standard signature functions [20] are normalized somewhat 
differently. The pomeron and reggeon trajectories are parametrized as a.p(t) = 1. 088 + 0. 25t 
and = 0.5 + 0.9t 

The differential cross section can be obtained from the corresponding amplitude as: 

= 77^|AWp°| 2 • (2-16) 
at 16ns 2 

The total cross section a can be obtained by integrating (12.161) over t 

rtmax(s) d& - 

a 77 ^ p o p o = —j dt , (2-17) 

where t m i n and t max are s-dependent kinematical limitations of t. 

In Figj3] we present the corresponding t-integrated cross section together with existing 
experimental data taken from 0,113] • Here the solid line includes the off-shell form factor 
while the dashed line does not. The vanishing of the VDM- Regge cross section at W 17 = 2m p 
is due to t m i n , t max limitations. It is obvious from Figj3]that the VDM- Regge model cannot 
explain the huge close-to-threshold enhancement. Up to now the origin of this enhancement 
remains unclear. The VDM-Regge model nicely describes the experimental data for W 17 > 
2.5 GeV. 



III. RESULTS 

The main ingredient for the calculation of photon fluxes is the charge form factor of the 
nucleus. In FigJDwe show the charge form factor of 197 Au calculated from the realistic charge 
density as measured in electron scattering off nuclei [3] • One can observe many oscillations 
characteristic for relatively sharp edge of the nucleus. This form factor is used to calculate 
flux of equivalent photons according to Eqs. fl2.2p and (12.31) . To illustrate the sensitivity of 
our calculation on details of the form factor, we also show a monopole form factor adjusted 
to the correct charge radius of the nucleus F ch ((f) = A 2 /(A 2 + q 2 ) , with A = 83MeV. 

In Figj5] we show the distribution of the cross section for the nucleus-nucleus scattering 
in photon-photon center-of-mass energy W 77 for both the low-energy component and high- 
energy VDM-Regge component. Below W 71 = 2 GeV the low-energy component dominates. 
The situation reverses above W 71 = 2 GeV. To study the high energy component one must 
impose an extra cut on W 71 , that is, the pp invariant mass M pp . However, at RHIC energies 
the high-energy tail of the photon spectrum is small, and the nuclear cross section drops 
quickly with increasing invariant mass of two-p mesons. By the solid lines we show the 
results for the realistic nuclear formfactor, while the dashed lines refer to the monopole form 
factor. The results obtained with the different form factors start to diverge at W 71 > 4 GeV. 
This is due to the increasing longitudinal momentum transfer involved for higher invariant 
mass of the produced system. The 7-factors at RHIC are not too large, so that the form- 
factor suppression is felt at such moderate invariant masses. The result for the point-like 
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FIG. 4: The modulus of the charge form factor F em (q) of the Au nucleus for realistic charge 
distribution (solid). For comparison we show the monopole form factor often used in practical 
applications (dashed). 



nucleus, which we show by the dotted line, overestimates the realistic cross section by more 
than an order of magnitude. 




FIG. 5: The Au + Au — ► Au + Au + p°p° cross section as a function of W 11 = M pp for the RHIC 
energy ^/s NN = 200 GeV. 

In FigJS]we show the distribution in impact parameter b = \bi— b 2 \ (see also FigJ2]). Again, 
we show distributions for the low- and high-energy components separately. Also shown are 



the distributions for the monopole form factor and for realistic charge density from [14 
One can clearly see different results for different approaches to calculate flux factors of 
equivalent photons. The sharp cutoff at b = 2Ra ~ 14 fm is precisely the ^-function coming 
from absorptive corrections. 

Finally in Figf7]we show distribution in rapidity-like variable Y. As far as small invariant 
masses of p°p° pairs dominate, we may refer to it as a rapidity of the p°p° pair. Compared 



S 




to the point-like case, the distribution obtained with realistic charge density is concentrated 
at midrapidities, and configurations when both p°'s are in very forward or both p°'s are in 
very backward directions are strongly damped compared to the case with point-like nucleus 
charges. At larger rapidities, one can again see a substantial difference between results 
obtained with an approximate monopole form factor and with the exact one calculated from 
realistic charge density. 




FIG. 7: The Au + Au — > Au + Au + p p cross section as a function of the rapidity of the p p pair 
Y for y/s NN = 200 GeV. In the left panel we show the decomposition into low- and high-energy 
contributions. The right panel shows the sensitivity to the nuclear form factor. 
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IV. CONCLUSIONS 



We have calculated, for the first time, realistic cross sections for exclusive p°p° produc- 
tion in ultrarelativistic heavy-ion collisions at RHIC in the framework of equivalent photon 
approximation. We have discussed uncertainties related to the way how the nuclear photon 
flux is calculated. We have used realistic charge densities to calculate the nuclear charge 
form factors. The absorption effects have been included. 

The low-energy part of the elementary 77 — > p°p° process has been parametrized and the 
parameters have been fitted to the e + e~ data while the high-energy part has been modeled 
in the vector-dominance Regge type model with parameters which are used to describe other 
hadronic processes. The model turned out to be consistent with the highest-energy data 
points (W ~ 3-4 GeV) from e + e~ collisions. 

It was shown that a realistic calculation of both ingredients is necessary to make reliable 
estimates of the nucleus-nucleus exclusive production of the p°p° pairs. Large cross sections, 
of the order of fraction of milibarn, have been found. The bulk of the cross section is, 
however, concentrated in low photon-photon energies (low p°p° invariant masses). Making 
cuts on higher p°p° invariant masses one can easily select high-energy component. The p° 
mesons, decaying into 7r + 7r~, can be measured e.g. by the STAR detector at RHIC. A Monte 
Carlo study is necessary to analyse feasibility of a measurement of the process discussed here. 

In the present analysis we have concentrated on processes with final nuclei in the ground 
state. It is very difficult, if not impossible, to measure such very forward/backward nuclei. 
The multiple Coulomb excitations associated with p°p° production may cause additional 
excitation of one or even both nuclei to the giant resonance region. The neutron emission 
from the giant resonances can be used then to tag the processes. The zero-degree calorimeters 
(ZDC) at RHIC can be used to measure such neutrons. We plan a detailed study of these 
processes in the future. 
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